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EDITORS'   PREFACE. 


THE  volume  called  Higher  Mathematics,  the  first  edition 
of  which  was  published  in  1896,  contained  eleven  Chapters  by 
eleven  authors,  each  chapter  being  independent  of  the  others, 
but  all  supposing  the  reader  to  have  at  least  a  mathematical 
training  equivalent  to  that  given  in  classical  and  engineering 
colleges.  ,The  publication  of  that  volume  is  now  discontinued 
and  the  chapters  are  issued  in  separate  form.  In  these  reissues 
it  will  generally  be  found  that  the  monographs  are  enlarged 
by  additional  articles  or  appendices  which  either  amplify  the 
former  presentation  or  record  recent  advances.  This  plan  of 
publication  has  been  arranged  in  order  to  meet  the  demand  of 
teachers  and  the  convenience  of  classes,  but  it  is  also  thought 
that  it  may  prove  advantageous  to  readers  in  special  lines  of 
mathematical  literature. 

It  is  the  intention  of  the  publishers  and  editors  to  add  other 
monographs  to  the  series  from  time  to  time,  if  the  call  for  the 
same  seems  to  warrant  it.  Among  the  topics  which  are  under 
consideration  are  those  of  elliptic  functions,  the  theory  of  num- 
bers, the  group  theory,  the.  calculus  of  variations,  and  non- 
Euclidean  geometry;  possibly  also  monographs  on  branches  of 
astronomy,  mechanics,  and  mathematical  physics  may  be  included. 
It  is  the  hope  of  the  editors  that  this  form  of  publication  may 
tend  to  promote  mathematical  study  and  research  over  a  wider 
field  than  that  which  the  former  volume  has  occupied. 

December,  1905.  ... 


AUTHOR'S   PREFACE. 


IN  republishing  this  short  treatise  in  book  form  the  author 
solicits  criticism  but  offers  no  apology.  The  type  of  the  book 
he  has  sought  to  imitate  is  that  shown  in  the  "  mathematical 
tracts  "  of  the  late  Sir  George  B.  Airy.  The  brevity  and  the  con- 
crete illustrations  of  these  "  tracts  "  have  served  very  effectively 
in  introducing  students  to  a  number  of  the  more  difficult  fields 
of  applied  mathematics;  and  it  is  hoped  that  this  treatise  will 
serve  a  similar  end. 

The  theory  of  probability  and  the  theory  of  errors  now  con- 
stitute a  formidable  body  of  knowledge  of  great  mathematical 
interest  and  of  great  practical  importance.  Though  developed 
largely  through  applications  to  the  more  precise  sciences  of  as- 
tronomy, geodesy,  and  physics,  their  range  of  applicability  extends 
to  all  of  the  sciences;  and  they  are  plainly  destined  to  play  an 
increasingly  important  r6le  in  the  development  and  in  the  appli- 
cations of  the  sciences  of  the  future.  Hence  their  study  is  not 
only  a  commendable  element  in  a  liberal  education,  but  some 
knowledge  of  them  is  essential  to  a  correct  understanding  of 
daily  events. 

No  special  novelty  of  presentation  is  claimed  for  this  work; 
but  the  reader  may  find  it  advantageous  to  know  that  a  definite 
plan  has  been  followed.  This  plan  consists  in  presenting  each 
principle,  first,  by  means  of  a  simple,  concrete  example;  passing, 
secondly,  to  a  general  statement  by  means  of  a  formula;  and, 
thirdly,  illustrating  applications  of  the  formula  by  concrete 
examples.  Great  pains  have  been  taken  also  to  secure  clear  and 
correct  statements  of  fundamental,  facts.  If  these  latter  are 
duly  understood,  the  student  needs  little  additional  aid;  if 
they  are  not  duly  understood,  no  amount  of  aid  will  forward  him. 

The  passage  from  the  elementary  concrete  to  the  advanced 
abstract  may  appear  to  be  abrupt  to  the  reader  in  some  cases. 
It  is  hoped,  however,  that  any  large  gaps  may  be  easily  bridged 
and  that  any  serious  difficulties  may  be  easily  overcome  by  means 
of  the  references  given  to  the  literature  of  the  subject.  In  any 
event  the  student  will  find  that  in  this,  as  in  all  of  the  more  ardu- 
ous sciences,  his  greatest  pleasure  and  his  highest  discipline  will 
come  from  bridging  such  gaps  and  from  surmounting  such  diffi- 
culties. 

WASHINGTON,  D.  C,  December,  1905. 
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PROBABILITY  AND  THEORY  OF  ERRORS, 


ART.  1.    INTRODUCTION. 

IT  is  a  curious  circumstance  that  a  science  so  profoundly 
mathematical  as  the  theory  of  probability  should  have  origi- 
nated in  the  games  of  chance  which  occupy  the  thoughtless  and 
the  profligate.*  That  such  is  the  case  is  sufficiently  attested 
by  the  fact  that  much  of  the  terminology  of  the  science  and 
many  of  its  familiar  illustrations  are  drawn  directly  from  the 
vocabulary  and  the  paraphernalia  of  the  gambler  and  the  trick- 
ster. It  is  somewhat  surprising,  also,  considering  the  antiquity 
of  games  of  chance,  that  formal  reasoning  on  the  simpler 
questions  in  probability  did  not  begin  before  the  time  of  Pascal 
and  Fermat.  Pascal  was  led  to  consider  the  subject  during  the 
year  1654  through  a  problem  proposed  to  him  by  the  Chevalier 
de  Mere,  a  reputed  gambler.f  The  problem  in  question  is 
known  as  the  problem  of  points  and  may  be  stated  as  follows : 
two  players  need  each  a  given  number  of  points  to  win  at  a 
certain  stage  of  their  game ;  if  they  stop  at  this  stage,  how  should 
the  stakes  be  divided?  Pascal  corresponded  with  his  friend 
Fermat  on  this  question  ;  and  it  appears  that  the  letters  which 
passed  between  them  contained  the  earliest  distinct  formulation 
of  principles  falling  within  the  theory  of  probability.  These 

*  The  historical  facts  referred  to  in  this  article  are  drawn  mostly  from  Tod- 
hunter's  History  of  the  Mathematical  Theory  of  Probability  from  the  time  of 
Pascal  to  that  of  Laplace  (Cambridge  and  London,  1865). 

f  "  Un  probleme  relatif  aux  jeux  de  hasard,  propose  aun  austere  janseniste 
par  un  homme  du  monde,  a  ete  1'origine  du  calcul  des  probabilites."  Poisson, 
Recherches  sur  la  Probabilite  des  Jugements  (Paris,  1837). 
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acute  thinkers,  however,  accomplished  little  more  than  a  correct 
start  in  the  science.  Each  seemed  to  rest  content  at  the  time 
with  the  approbation  of  the  other.  Pascal  soon  renounced 
such  mundane  studies  altogether ;  Fermat  had  only  the  scant 
leisure  of  a  life  busy  with  affairs  to  devote  to  mathematics; 
and  both  died  soon  after  the  epoch  in  question, — Pascal  in 
1662,  and  Fermat  in  1665. 

A  subject  which  had  attracted  the  attention  of  such  dis- 
tinguished mathematicians  could  not  fail  to  excite  the  interest 
of  their  contemporaries  and  successors.  Amongst  the  former 
Huygens  is  the  most  noted.  He  has  the  honor,  of  publishing 
the  first  treatise*  on  the  subject.  It  contains  only  fourteen 
propositions  and  is  devoted  entirely  to  games  of  chance,  but  it 
gave  the  best  account  of  the  theory  down  to  the  beginning  of 
the  eighteenth  century,  when  it  was  superseded  by  the  more  elab- 
orate works  of  James  Bernoulli,f  Montmort^  and  De  Moivre.g 
Through  the  labors  of  the  latter  authors  the  mathematical 
theory  of  probability  was  greatly  extended.  They  attacked, 
quite  successfully  in  the  main,  the  most  difficult  problems ; 
and  great  credit  is  due  them  for  the  energy  and  ability  dis- 
played in  developing  a  science  which  seemed  at  the  time  to 
have  no  higher  aim  than  intellectual  di version. |  Their  names, 
undoubtedly,  with  one  exception,  that  of  Laplace,  are  the  most 
important  in  the  history  of  probability. 

Since  the  beginning  of  the  eighteenth  century  almost  every 
mathematician  of  note  has  been  a  contributor  to  or  an  expos- 
itor of  the  theory  of  probability.  Nicolas,  Daniel,  and  John 
Bernoulli,  Simpson,  Euler,  d'Alembert,  Bayes,  Lagrange,  Lam- 
bert, Condorcet,  and  Laplace  are  the  principal  names  which 
figure  in  the  history  of  the  subject  during  the  hundred  years 

*De  Ratiociniis  in  Ludo  Aleae,  1657. 

f  Ars  Conjectandi,  1713. 

JEssai  d'Analyse  sur  les  Jeux  de  Hazards,  1708. 

§The  Doctrine  of  Chances,  1718. 

[Todhunter  says  of  *Montmort,  for  example,  "  In  1708  he  published  his 
work  on  Chances,  where  with  the  courage  of  Columbus  he  revealed  anew  world 
to  Mathematicians." 
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ending  with  the  first  quarter  of  the  nineteenth  century.  Of 
contributions  from  this  brilliant  array  of  mathematical  talent, 
the  Th£orie  Analytique  des  Probability  of  Laplace  is  by  far 
the  most  profound  and  comprehensive.  It  is,  like  his  Me- 
canique  Celeste  in  dynamical  astronomy,  still  the  most  elabo- 
rate treatise  on  the  subject.  An  idea  of  the  grand  scale  of  the 
work  in  its  present  form*  may  be  gained  by  the  facts  that  the 
non-mathematical  introduction-)-  covers  about  one  hundred  and 
fifty  quarto  pages ;  and  that,  in  spite  of  the  extraordinary 
brevity  of  mathematical  language,  the  pure  theory  and  its  ac- 
cessories and  applications  require  about  six  hundred  and  fifty 
pages. 

From  the  epoch  of  Laplace  down  to  the  present  time  the 
extensions  of  the  science  have  been  most  noteworthy  in  the 
fields  of  practical  applications,  as  in  the  adjustment  of  obser- 
vations, and  in  problems  of  insurance,  statistics,  etc.  Amongst 
the  most  important  of  the  pioneers  in  these  fields  should 
be  mentioned  Poisson,  Gauss,  Bessel,  and  De  Morgan.  Nu- 
merous authors,  also,  have  done  much  to  simplify  one  or  an- 
other branch  of  the  subject  and  thus  bring  it  within  the  range 
of  elementary  presentation.  The  fundamental  principles  of 
the  theory  are,  indeed,  now  accessible  in  the  best  text-books 
on.algebra :  and  there  are  many  excellent  treatises  on  the  pure 
theory  and  its  various  applications. 

Of  all  the  applications  of  the  doctrine  of  probability  none 
is  of  greater  utility  than  the  theory  of  errors.  In  astronomy, 
geodesy,  physics,  and  chemistry,  as  in  every  science  which  at- 
tains precision  in  measuring,  weighing,  and  computing,  a 
knowledge  of  the  theory  of  errors  is  indispensable.  By  the  aid 
of  this  theory  the  exact  sciences  have  made  great  progress  dur- 

*The  form  of  the  third  edition  published  in  1820,  and  of  Vol.  VII  of  the 
complete  works  of  Laplace  recently  republished  under  the  auspices  of  the 
Academic  des  Sciences  by  Gauthier-Villars.  This  Vol.  VII  bears  the  date  1886. 

f  "  Cette  Introduction,"  writes  Laplace,  "est  le  developpement  d'une  Le9on 
sur  les  Probabilites,  que  je  donnai  en  1795,  aux  Ecoles  Normales,  ou  je  fus  ap- 
pele  comme  professeur  de  Mathematiques  avec  Lagrange,  par  un  decrct  de  la 
Convention  nationale." 
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ing  the  nineteenth  century,  not  only  in  the  actual  determination 
of  the  constants  of  nature,  but  also  in  the  fixation  of  clear 
ideas  as  to  the  possibilities  of  future  conquests  in. the  same  di- 
rection. Nothing,  for  example,  is  more  satisfactory  and  in- 
structive in  the  history  of  science  than  the  success  with  which 
the  unique  method  of  least  squares  has  been  applied  to  the 
problems  presented  by  the  earth  and  the  other  members  of  the 
solar  system.  So  great,  in  fact,  are  the  practical  value  and 
theoretical  importance  of  the  method  of  least  squares,  that  it  is 
frequently  mistaken  for  the  whole  theory  of  errors,  and  is 
sometimes  regarded  as  embodying  the  major  part  of  the  doc- 
trine of  probability  itself. 

As  may  be  inferred  from  this  brief  sketch,  the  theory  of 
probability  and  its  more  important  applications  now  constitute 
an  extensive  body  of  mathematical  principles  and  precepts. 
Obviously,  therefore,  it  will  be  impossible  within  the  limits  of 
a  single  condensed  monograph  to  do  more  than  give  an  out- 
line of  the  salient  features  of  the  subject.  It  is  hoped,  how- 
ever, in  accordance  with  the  general  plan  of  the  volume,  that 
such  outline  will  prove  suggestive  and  helpful  to  those  who 
may  come  to  the  science  for  the  first  time,  and  also  to  those 
who,  while  somewhat  familiar  with  the  difficulties  to  be  over- 
tome,  have  not  acquired  a  working  knowledge  of  the  subject. 
Effort  has  been  made  especially  to  clear  up  the  difficulties  of 
the  theory  of  errors  by  presenting  a  somewhat  broader  view  of 
the  elements  of  the  subject  than  is  found  in  the  standard 
treatises,  which  confine  attention  almost  exclusively  to  the 
method  of  least  squares.  This  chapter  stops  short  of  that 
method,  and  seeks  to  supply  those  phases  of  the  theory  which 
are  either  notably  lacking  or  notably  erroneous  in  works 
hitherto  published.  It  is  believed,  also,  that  the  elements  here 
presented  are  essential  to  an  adequate  understanding  of  the 
well-worked  domain  of  least  squares.* 

*The  author  has  given  a  brief  but  comprehensive  statement  of  the  method 
of  least  squares  in  the  volume  of  Geographical  Tables  published  by  the  Smith- 
sonion  Institution,  1894. 
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The  formulas  and  results  of  the  theory  of  permutations 
and  combinations  are  often  needed  for  the  statement  and  so- 
lution of  problems  in  probabilities.  This  theory  is  now  to  be 
found  in  most  works  on  algebra,  and  it  will  therefore  suffice 
here  to  state  the  principal  formulas  and  illustrate  their  mean- 
ing by  a  few  numerical  examples. 

The  number  of  permutations  of  n  things  taken  r  in  a  group 
is  expressed  by  the  formula 

(n)r  =  n(n  -  i)(«  -  2)  .  .  .  (n  -  r  +  i).  (i) 

Thus,  to  illustrate,  the  number  of  ways  the  four  letters  a,  6, 

c,  d  can  be  arranged  in  groups  of  two  is  4 .  3  =  12,  and  the  groups 
are 

abt     ba,     ac,     ca,     ad,     da,     be,     cb,     bd,     db,     cd,     dc. 
Similarly,  the  formula  gives  for 
n  =    3  and  r  —  2,         (3),  =  3.2  =6, 

n  =    7    "     r  =  3,         (7)3  =  7.6.5  =210, 

n  =  10    "     r  =  6,       (10)6=10.9.8.7.6.5  =  151200. 
The  results  which  follow  from  equation  (i)  when  n  and  r 
do  not  exceed  10  each  are  embodied  in  the  following  table : 
VALUES  OF  PERMUTATIONS. 


10 

9 

8 

7 

6 

5 

4 

3 

2 

z 

T 

10 

9 

8 

7 

6 

5 

4 

3 

2 

I 

2 

90 

72 

c6 

42 

30 

20 

12 

6 

2 

3 

720 

504 

336 

210 

120 

60 

24 

6 

4 

5040 

3024 

1680 

840 

36O 

1  20 

24 

5 

30240 

15120 

6720 

2520 

720 

120 

6 

i  5  i  200 

60480 

20160 

5040 

720 

7 

604800 

181440 

40320 

5040 

8 

1814400 

362880 

40320 

9 

3628800 

362880 

10 

3628800 

Sp 

9864100 

986409 

109600 

13699 

1956 

325 

64 

15 

4 

I 

The  use  of  this  table  is  obvious.  Thus,  the  number  of  per- 
mutations of  eight  things  in  groups  of  five  each  is  found  in  the 
fifth  line  of  the  column  headed  with  the  number  8.  It  will  be 
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noticed  that  the  last  two  numbers  in  each  column  (excepting 
that  headed  with  i)  are  the  same.  This  accords  with  the  for- 
mula,  which  gives  for  the  number  of  permutations  of  n  things 
in  groups  of  n  the  same  value  as  for  n  things  in  groups  of 
(;/  —  i).  It  will  also  be  remarked  that  the  last  number  in  each 
column  of  the  table  is  the  factorial,  n  !,  of  the  number  n  at  the 
head  of  the  column.  For  example,  in  the  column  under  7,  the 
last  number  is  5040  =  1.2.3.4.5.6.7  =  7!. 

The  total  number  of  permutations  of  n  things  taken  singly, 
in  groups  of  two,  three,  etc.,  is  found  by  summing  the  numbers 
given  by  equation  (i)  for  all  values  of  r  from  I  jto  n.  Calling 
this  total  or  sum  Sp,  it  will  be  given  by 

S,  =  2(n\.  (2) 

To  illustrate,  suppose  n  =  3,  and,  to  fix  the  ideas,  let  the 
three  things  be  the  three  digits  I,  2,  3.  Then  from  the  above 
table  it  is  seen  that  Sf  =  3  +  6  -J-  6  —  15  ;  or,  that  the  number 
of  numbers  (all  different)  which  can  be  formed  from  those  dig- 
its is  fifteen.  These  numbers  are  I,  2,  3  ;  12,  13,  21,  23,  31,  32  ; 
123,  132,  213,  231,  312,  321. 

The  values  of  Sj,  for  n  =  I,  2,  .  .  .  10  are  given  under  the 
corresponding  columns  of  the  above  table.  But  when  n  is 
large  the  direct  summation  indicated  by  (2)  is  tedious,  if  not 
impracticable.  Hence  a  more  convenient  formula  is  desirable. 
To  get  this,  observe  that  (i)  may  be  written 


if  r  is  restricted  to  integer  values  between  I  and  (n  —  i),  both 
inclusive.  This  suffices  to  give  all  terms  which  appear  in  the 
right-hand  member  of  (2),  since  the  number  of  permutations 
for  r  =  (n  —  i)  is  the  same  as  for  r  =  n.  Hence  it  appears 
that 

"!        nl  nl 


COMBINATIONS.  13 

But  as  n  increases,  the  series  by  which  n !  is  here  multiplied 
approximates  rapidly  towards  the  base  of  natural  logarithms ; 
that  is,  towards 

^  =  2.7182818+,         log  e  =  0.4342945. 
Hence  for  large  values  of  n 

Sp  =  n\  e,  approximately.*  (3) 

To  get  an  idea  of  the  degree  of  approximation  of  (3),  sup- 
pose n  =  9.  Then  the  computation  runs  thus  (see  values  in 

the  above  table) : 

log 

9!  =  362880         5. 5  597630 
e  0.4342945 

gle  =  986410         5-9940575 
Sf  =  986409  by  equation  (2). 

The  error  in  this  case  is  thus  seen  to  be  only  one  unit,  or 
about  one-millionth  of  S^.f 

Prob.  i.  Tabulate  a  list  of  the  numbers  of  three  figures  each 
which  can  be  formed  from  the  first  five  digits  i,  .  .  .  5.  How  many 
numbers  can  be  formed  from  the  nine  digits  ? 

Prob.  2.  Is  Sj,  always  an  odd  number  for  n  odd  ?  Observe 
values  of  Sp  in  the  table  above. 

ART.  3.    COMBINATIONS. 

In  permutations  attention  is  given  to  the  order  of  arrange- 
ment of  the  things  considered.  In  combinations  no  regard  is 
paid  to  the  order  of  arrangement.  Thus,  the  permutations  of 
the  letters  a,  b,  c,  d  in  groups  of  three  are 

(abc)  (abd)  bac  bad  acb  (acd]  cab  cad 
adb  adc  dab  dac  bca  (bed)  cba  cbd 
bda  bdc  dba  dbc  cda  cdb  dca  deb 

*  See  Art.  6  for  a  formula  for  computing '«!  when  «  is  a  large  number. 

f  When  large  numbers  are  to  be  dealt  with,  equations  (i)'  and  (3)  are  easily 
managed  by  logarithms,  especially  if  a  table  of  values  of  log  («!)  is  available. 
Such  tables  are  given  to  six  places  in  De  Morgan'?  treatise  on  Probability  in 
the  Encyclopaedia  Metropolitana,  and  to  five  places  in  Shortrede's  Tables 
(Vol.  I,  1849). 
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But  if  the  order  of  arrangement  is  ignored  all  of  these  are 
seen  to  be  repetitions  of  the  groups  enclosed  in  parentheses, 
namely,  (abc\  (abd),  (acd),  (bcd\  Hence  in  this  case  out  of 
twenty-four  permutations  there  are  only  four  combinations. 

A  general  formula  for  computing  the  number  of  combina- 
tions of  n  things  taken  in  groups  of  r  things  is  easily  derived. 
For  the  number  of  permutations  of  n  things  in  groups  of  r  is 
by  (i)  of  Art.  2 

(n)r  =  n(n  -  i)(»  -  2)  ...(»-  r  +  i)  ; 

and  since  each  group  of  r  things  gives  I  .  2  .  3  .  .  .  r  =  r\  per- 
mutations, the  number  of  combinations  must  be  the  quotient 
of  (n)r  by  r\.  Denote  this  number  by  C(n)r.  Then  the  gen- 
eral formula  is 


This  formula  gives,  for  example,  in  the  case  of  the  four  let- 
ters a,  by  ct  d  taken  in  groups  of  three,  as  considered  above, 


Multiply  both  numerator  and  denominator  of  the  right-hand 
member  of  (i)  by  (n  —  r)l     The  result  is 


(I)' 


which  shows  that  the  number  of  combinations  of  n  things  in 
groups  of  r  is  the  same  as  the  number  of  combinations  of  n 
things  in  groups  of  (n  —  r).  Thus,  the  number  of  combina- 
tions of  the  first  ten  letters  a,  b,  c  .  .  .j  in  groups  of  three  or 
seven  is 

10  ! 


The  following  table  gives  the  values  C(n)r  for  all  values  of 
n  and  r  from  I  to  10, 

The  mode  of  using  this  table  is  evident.  For  example,  the 
number  of  combinations  of  eight  things  in  sets  of  five  each  is 
.found  on  the  fifth  line  of  the  column  headed  8  to  be  56. 
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VALUES  OF  COMBINATIONS. 


15 


10 

9 

8 

7 

6 

5 

4 

3 

2 

I 

I 

10 

9 

8 

7 

6 

s 

4 

3 

2 

i 

2 

45 

36 

28 

21 

IS 

10 

6 

1 

I 

3 

120 

84 

56 

35 

20 

10 

4 

i 

4 

2IO 

126 

70 

35 

IS 

5 

I 

1 

252 
2IO 

126 

84 

56 
28 

21 

7 

6 

i 

i 

I 

120 

45 

36 
9 

8 

i 

i 

9 

10 

i 

10 

i 

St 

1023 

5ii 

255 

127 

63 

3i 

15 

7 

3 

It  will  be  observed  that  the  numbers  in  any  column  show 
a  maximum  value  when  n  is  even  and  two  equal  maximum 
values  when  n  is  odd.  That  this  should  be  so  is  easily  seen 
from  (i)',  which  shows  that  C(n)r  will  be  a  maximum  for  any 
value  of  n  when  r\  (n  —  r)  !  is  a  minimum.  For  n  even  this  is 
a  minimum  for  r  =  ^n  ;  while  for  n  odd  it  has  equal  minimum 
values  for  r  =  %(n  —  i)  and  r  —  %(n  +  i).  Thus, 


maximum  of  C(n\  = 


for  n  even, 


(2) 


for  »  odd. 


The  total  number  of  combinations  of  n  things  taken  singly, 
in  groups  of  two,  three,  etc.,  is  found  by  summing  the  numbers 
given  by  (i)  for  all  values  of  r  from  I  to  n  both  inclusive. 
Calling  this  total  or  sum  Sc, 

Sc  =  2C(n)r. 

The  same  sum  will  also  come  from  (i)'by  giving  to  rail  values 
from  i  to  (n  —  i),  both  inclusive,  summing  the  results,  and  in- 
creasing their  aggregate  by  unity.     Thus  by  either  process 
,    n(n  -  i)   ,    n(n  -  i](n  -3) 


I  .  2 


1.2.3 
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The  second  member   of   this  equation  is  evidently  equal   to 
(i  +  0"  T  *•     Hence 

Se=3C(n),  =  2*-i.  (3) 

The  values  of  Sc  for  values  of  n  and  r  from  I  to  10  are  given 
under  the  corresponding  columns  of  the  above  table. 

Prob.  3.  How  many  different  squads  of  ten  men  each  can  be 
formed  from  a  company  of  100  men  ? 

Prob.  4  How  many  triangles  are  formed  by  six  straight  lines 
each  of  which  intersects  the  other  five  ? 

Prob.  5.  Examine  this  statement  :  "  In  dealing  a  pack  of  cards 
the  number  of  hands,  of  thirteen  cards  each,  which^can  be  produced 
is  635  013  559  600.  But  in  whist  four  hands  are  simultaneously  held,. 
and  the  number  of  distinct  deals  .  .  .  would  require  twenty-eight 
figures  to  express  it."  * 

Prob.  6.  Assuming  combination  always  possible,  and  disregarding 
the  question  of  proportions,  find  how  many  different  substances 
could  be  produced  by  combining  the  seventy-three  chemical  ele- 
ments. 

ART.  4.     DIRECT  PROBABILITIES. 

If  it  is  known  that  one  of  two  events  must  occur  in  any 
trial  or  instance,  and  that  the  first  can  occur  in  a  ways  and  the 
second  in  b  ways,  all  of  which  ways  are  equally  likely  to  hap- 
pen, then  the  probability  that  the  first  will  happen  is  expressed 
mathematically  by  the  fraction  a/(a-\-b),  while  the  probability 
that  the  second  will  happen  is  b/(a  -f-  b).  Such  events  are  said 
to  be  mutually  exclusive.  Denote  their  probabilities  by  p  and 
q  respectively.  Then  there  result 


the  last  equation  following  from  the  first  two  and  being  the 
mathematical  expression  for  the  certainty  that  one  of  the  two 
events  must  happen. 

Thus,  to  illustrate,  in  tossing  a  coin  it  must  give  "  head  "  or 
"  tail"  ;  a  =  b  =  I,  and  p  =  q  =  1/2.  Again,  if  an  urn  contain 
a  =  5  white  and  b  =  8  black  balls,  the  probability  of  drawing 

*  Jevons,  Principles  of  Science,  New  York,  1874,  p.  217. 
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a  white  ball  in  one  trial  is  /  =  5/13  and  that  of  drawing  a 
black  one  q  =  8/13. 

Similarly,  if  there  are  several  mutually  exclusive  events 
which  can  occur  in  a,  b,  c  .  .  .  ways  respectively,  their  probabil- 
ities /,  q,  r  .  .  .  are  given  by 

_  a  —  ^  _  c 

+  .  .  .'  T  =  a  -\-b-\-  *+...' 

(2) 

..=  l, 

For  example,  if  an  urn  contain  a  =  4  white,  #=5  black, 
and  c  =  6  red  balls,  the  probabilities  of  drawing  a  white,  black, 
and  red  ball  at  a  single  trial  are  ^  =  4/15,  ^=5/15,  and 
r  =  6/1$,  respectively. 

Formulas  (i)  and  (2)  may  be  applied  to  a  wide  variety  of 
cases,  but  it  must  suffice  here  to  give  only  a  few  such.  As  a 
first  illustration,  consider  the  probability  of  drawing  at  random 
a  number  of  three  figures  from  the  entire  list  of  numbers  which 
can  be  formed  from  the  first  seven  digits.  A  glance  at  the 
table  of  Art.  I  shows  that  the  symbols  of  formula  (i)  have  in 
this  case  the  values  #—210,  and  #+£=  13699.  Hence 
b  =  13489,  and  /  =  210/13699  ;  that  is,  the  probability  in  ques- 
tion is  about  1/65. 

Secondly,  what  is  the  probability  of  holding  in  a  hand  of 
whist  all  the  cards  of  one  suit  ?  Formula  (i)  of  Art.  3  shows 
that  the  number  of  different  hands  of  thirteen  cards  each  which 
may  be  formed  from  a  pack  of  fifty-two  cards  is 


and  the  probability  required  is  the  reciprocal  of  this  number.. 
The  probability  against  this  event  is,  therefore,  very  nearly 
unity. 

Thirdly,  consider  the  probabilities  presented  by  the  case  of 
an  urn  containing  4  white,  5  black,  and  6  red  balls,  from  which 
at  a  single  trial  three  balls  are  to  be  drawn.  Evidently  the 
triad  of  balls  drawn  may  be  all  white,  all  black,  all  red,  partly 
white  and  black,  partly  white  and  red,  partly  black  and  red,  or 


—  (io-f-2o)  =  i35  =f 
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•one  each  of  the  white,  black,  and  red.  There  are  thus  seven 
-different  probabilities  to  be  taken  into  account.  The  theory 
of  combinations  shows  (see  equation  (i),  Art,  3)  that  the  total 
number  of 

White  triads 

Black  triads 

Red  triads 

White  and  black  triads 

White  and  red  triads 

Elack  and  red  triads 

White,  black,  and  red  triads  =  4.5.6  =  120  =  g 

Sum  =  455 

The  total  number  of  these  triads  is  455,  and  is,  as  it  should 
;be,  the  number  of  combinations  in  groups  of  three  each  of  the 
-whole  number  of  balls.     Hence  formulas  (2)  give  the  seven 
different  probabilities   which   follow,   using  the    initial  letters 
-zv,  b,  r  to  indicate  the  colors  represented  in  a  triad  : 
For  a  triad  www  p  =      4/455, 

"    "     "      bbb  q  =    10/455, 

"     "     "      rrr  r  =    20/455, 

«    «     u      Wwb  or  wbb  s=    70/455, 

"    "     "      wwr  or  wrr  /  ==    96/455, 

"    "     "      bbr  or  brr  u  =  135/455, 

«    «     .<      wfrr  v  —  120/455. 

Prob.  7.  When  three  dice  are  thrown  together,  what  is  the  prob- 
ability that  the  throw  will  be  greater  than  9  ? 

Prob.  8.  Write  down  a  literal  formula  for  the  probabilities  of  the 
several  possible  triads  considered  in  the  above  question  of  the  balls, 
supposing  the  numbers  of  white,  black,  and  red  balls  to  be  /,  m,  n, 
.  respectively. 
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ART.  5.    PROBABILITY  OF  CONCURRENT  EVENTS. 

If  the  probabilities  of  two  independent  events  are/,  and 
/2,  respectively,  the  probability  of  their  concurrence  in  any 
single  instance  is  /,/>,.  Thus,  suppose  there  are  two  urns 
£/i  and  £7a,  the  first  of  which  contains  al  white  and  £,  black 
balls,  and  the  second  <za  white  and  bt  black  balls.  Then  the 
probability  of  drawing  a  white  ball  from  £7,  is/,  =  al/(al  +  b^ 
while  that  of  drawing  a  white  ball  fromxC/"2  is/a  =  <za/(tfa+  £,)- 
The  total  number  of  different  pairs  of  balls  which  can  be  formed 
from  the  entire  number  of  balls  is  (<z,  -f-  b^(a%  +  £,)•  Of  these 
pairs  a^a%  are  favorable  to  the  concurrence  of  white  in  simul- 
taneous or  successive  drawings  from  the  two  urns.  Hence  the 
probability  of  a  concurrence  of 

white  with  white  =  alaz/(a1  -f  b^}(at  +  £,), 

white  with  black  =  (a&  +  a^bl)/(al  +  £,)(#,  +  &,), 

black  with  black  =  b,bz/(a,  +  b,)(az  -f  b^\ 

and  the  sum  of  these  is  unity,  as  required  by  equations  (2)  of 
Article  4. 

In  general,  if  /,,  /„,  pt  .  .  .  denote  the  probabilities  of  several 
independent  events,  and  P  denote  the  probability  of  theii 
concurrence, 

^  =  AAA--.  (i) 

To  illustrate  this  formula,  suppose  there    is  required  the 
probability  of  getting  three  aces  with  three  dice  thrown  simul- 
taneously.    In  this  case/,  —  /a  =/3  =  1/6  and 
P=(i/6)3=  1/216. 

Similarly,  if  two  dice  are  thrown  simultaneously  the  proba- 
bility that  the  sum  of  the  numbers  shown  will  be  1 1  is  2/36; 
and  the  probability  that  this  sum  1 1  will  appear  in  two  succes- 
sive throws  of  the  same  pair  of  dice  is  4/36.36. 

The  probability  that  the  alternatives  of  a  series  of  events, 
will  concur  is  evidently  given  bv 

Q  =  <?><?,<?, . . .  =(i  -  AX'  ~  A)0  -A)  •  •  -  (2) 

Thus,  in  the  case  of  the  three  dice  mentioned  above,  the 

probability  that  each  will  show  something  other  than  an  ace  is 
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#,  =  q%  =  qt  =  5/6,  and  the  probability  that  they  will  concur  in 
this  is  Q  =  125/216. 

Many  cases  of  interest  occur  for  the  application  of  (i)  and 
(2).  One  of  the  most  important  of  these  is  furnished  by  suc- 
cessive trials  of  the  same  event.  Consider,  for  example,  what 
may  happen  in  n  trials  of  an  event  for  which  the  probability 
is  p  and  agajnst  which  the  probability  is  q.  The  probability 
that  the  event  will  occur  every  time  is/n.  The  probability  that 
the  event  will  occur  (n  —  i)  times  in  succession  and  then  fail  is 
J)n  ~  lq.  But  if  the  order  of  occurrence  is  disregarded  this  last 
combination  may  arrive  in  n  different  ways  ;  so  that  the  prob- 
ability that  the  event  will  occur  (n  —  i)  times  and  fail  once  is 
npn~*q.  Similarly,  the  probability  that  the  event  will  happen 
{11  —  2)  times  and  fail  twice  is  \n(n  —  i)  pll~  y ;  etc-  That  is, 
the  probabilities  of  the  several  possible  occurrences  are  given  by 
the  corresponding  terms  in  the  development  of  (p  -j-  q)n. 

By  the  same  reasoning  used  to  get  equations  (2)  of  Art. 
3  it  may  be  shown  that  the  maximum  term  in  the  expansion 
of  (p-\-q)n  is  that  in  which  the  exponent  m,  say,  of  q  is 
the  whole  number  lying  between  (n-\-  i)q  —  I  and  (n  +  \)q. 
In  other  words,  the  most  probable  result  in  n  trials  is  the 
occurrence  of  the  event  (n  —  m)  times  and  its  failure  m 
times.  When  n  is  large  this  means  that  the  most  probable  of 
all  possible  results  is  that  in  which  the  event  occurs  n  —  nq 
=  n(\  —  q)  =  np  times  and  fails  nq  times.  Thus,  if  the  event 
be  that  of  throwing  an  ace  with  a  single  die  the  most  probable 
of  the  possible  results  in  600  throws  is  that  of  100  aces  and 
500  failures. 

Since  qn  is  the  probability  that  the  event  will  fail  every  time 
in  n  trials,  the  probability  that  it  will  occur  at  least  once  in  n 
trials  is  I  —  qn.  Calling  this  probability  r* 

r=  i  -?*=  i  -(i  -/)«.  (3) 

If  r  in  this  equation  be  replaced  by  1/2,  the  corresponding 
value  of  n  is  the  number  of  trials  essential  to  render  the 

*  See  Poisson's  Probabilite  des  Jugements,  pp.  40.  41. 
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chances  even  that  the  event  whose  probability  is/  will  occur 
at  least  once.     Thus,  in  this  case,  the  value  of  n  is  given  by 


log  (i-p] 

This  shows,  for  example,  if  the  event  be  the  throwing  of  double 
sixes  with  two  dice,  for  which  /  =  1/36,  that  the  chances  are 
even  (r  =  1/2)  that  in  25  throws  (n  —  24.614  by  the  formula) 
double  sixes  will  appear  at  least  once. 

Equation  (3)  shows  that  however  small/  may  be,  so  long  as 
it  is  finite,  n  may  be  taken  so  large  as  to  make  r  approach  in- 
definitely near  to  unity  ;  that  is,  n  may  be  so  large  as  to  render 
it  practically  certain  that  the  event  will  occur  at  least  once. 

When  n  is  large 

n(n  —  i)    .      n(n  —  i}(n—  2)  , 

-^TjV-      1.2.3   y  +  --» 


=  e~np  approximately. 
Thus  an  approximate  value  of  r  is 

r  =  i  -  e  -  *>,         \oge  =  O.4342495-  (4) 

This  formula  gives,  for  example,  for  the  probability  of  drawing 
the  ace  of  spades  from  a  pack  of  fifty-two  cards  at  least  once  in 
104  trials  r  =  i  —  e~*  =:  0.865,  while  the  exact  formula  (3) 
gives  0.867. 

Similarly,  the  probability  of  the  occurrence  of  the  event  at 
least  /  times  in  n  trials  will  be  given  by  the  sum  of  the  terms 
Of  (p  _|_  qY  from  pn  up  to  that  in  ptqn~t  inclusive.  This  proba- 
bility must  be  carefully  distinguished  from  the  probability  that 
the  event  will  occur  /  times  only  in  the  n  trials,  the  latter  being 
expressed  by  the  single  term  in  /V"'. 

Prob.  9.  Compare  the  probability  of  holding  exactly  four  aces  in 
five  hands  of  whist  with  the  probability  cf  holding  at  least  four  aces 
in  the  same  number  of  hands. 

Prob.  10.  What  is  the  probability  of  an  event  if  the  chances  are 
even  that  it  occurs  at  least  once  in  a  million  trials  ?  See  equation  (4). 
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ART.  6.    BERNOULLI'S  THEOREM. 

Denote  the  exponents  of/  and  q  in  the  maximum  term  of 
(/  +  !?)*  by  /*  anc*  m  respectively,  and  denote  this  term  by  7". 
Then 

n(n  —  \\(n  —  2)  .  .  .  (yw  -4-  i 

-- 


m  }Ji\  m 

As  shown  in  Art.  5,  /*  in  this  formula  is  the  greatest  whole 
number  in  («  +  i)/,  and  ?«  the  greatest  whole  number  in 
(n  -\-  \)q\  so  that  when  n  is  large,  /t  and  m  are  sensibly  equal 
to  np  and  #^  respectively. 

The  direct  calculation  of  T  by  (i)  is  impracticable  when  n 
is  large.  To  overcome  this  difficulty  the  following  expression 
is  used  :  * 

(2) 


log  e  =  0.4342495,     log  2n  =  0.7981799. 

This  expression  approaches  nne~n  ^2nn  as  a  limit  with  the 
increase  of  #,  and  in  this  approximate  form  is  known  as  Stir- 
ling's theorem.  Although  a  rude  approximation  to  n  \  for 
small  values  of  n  this  theorem  suffices  in  nearly  all  cases 
wherein  such  probabilities  as  T  are  desired.  Making  use  of 
the  theorem  in  (i)  it  becomes 


That  this  formula  affords  a  fair  approximation  even  when 
n  is  small  is  seen  from  the  case  of  a  die  thrown  12  times.  The 
probability  that  any  particular  face  will  appear  in  one  throw  is 
/  =  1/6,  whence  q  =  5/6;  and  the  most  probable  result  in  12 
throws  is  that  in  which  the  particular  face  appears  twice  and 
fails  to  appear  ten  times.  The  probability  of  this  result  com- 
puted from  (3)  is  0.309,  while  the  exact  formula  (i)  gives  0.296. 

The  probability  that  the  event  will  occur  a  number  of  times 

*  This  expression  is  due  to  Laplace,  Th6orie  Analytique  des  Probabilit6s. 
See  also  De  Morgan's  Calculus,  pp.  600-604. 
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comprised  between  (yu  —  a)  and  (yu  +  a)  in  n  trials  is  evidently 
expressed  by  the  sum  of  the  terms  in  (p  -f  q)n  for  which  the 
exponent  of  /  has  the  specified  range  of  values.  Calling  this 
probability  R,  putting 

yu  =  up  -\-  u,     and     m  =  nq  —  uy 

and  using  Stirling's  theorem  (which  implies  that  n  is  a  large 
number),* 


R  =  2  i  4- 

V27tnpq\         npl  nq 

very  nearly  ;  and  the  summation  is  with  respect  to  u  from 
u  =  —  a  to  u  =  -{-  a.  But  expansion  shows  that  the  natural 
logarithm  of  the  product  of  the  two  binomial  factors  in  this 
equation  is  approximately  —  u*/2npq.  Hence 


V2itnpq 

and,  since  n  is  supposed  large,  this  may  be  replaced  by  a  definite 
integral,  putting 


dz  =  i/V2npq,     and    £  =  u*/2npq. 
Thus 

+  a 

(4) 


+  a/  t/2^"  a/  Vw~q 

=^  fe-^dz  =  -2=  [r* 

VnJ  VnJ 


This  equation  expresses  the  theorem  of  James  Bernoulli, 
given  in  his  Ars  Conjectandi,  published  in  1713. 

The  value  of  the  right-hand  member  of  (4)  varies,  as  it 
should,  between  o  and  i,  and  approaches  the  latter  limit  rap- 
idly as  z  increases.  Thus,  writing  for  brevity 


/=-T= 


*  See  Bertrand,  Calcul  des  Probabilites,  Paris,  1889,  for  an  extended  discus- 
sion of  the  questions  considered  in  this  Article. 
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the  following  table  shows  the  march  of  the  integral 


z 

' 

z 

I 

* 

/ 

0.00 
•25 

0.000 

.276 

0-75 
1.  00 

0.711 

.843 

1.50 

1-75 

0.966 
.987 

•50 

.520 

1.25 

.923 

2.00 

•995 

To  illustrate  the  use  of  (4),  suppose  there  is  required  the 
probability  that  in  6000  throws  of  a  die  the  ace  will  appear  a 
number  of  times  which  shall  be  greater  than  1/6  X  6000  —  10 
and  less  than  1/6  X  6000  ~\-  10,  or  a  number  of  times  lying 
between  990  and  1010.  In  this  case  a  =  10,  n  =  6000,  /  =  1/6, 
^  =  5/6.  Thus,  a/V2npq  =  10/^2  .  6000  .  1/6  .  5/6  =  0.245. 
Hence,  by  (4)  and  the  table,  R  =  0.27. 

Prob.  ii.  If  the  ratio  of  males  to  females  at  birth  is  105  to  100, 
what  is  the  probability  that  in  the  next  10,000  births  the  number  of 
males  will  fall  within  two  per  cent  of  the  most  probable  number? 

Prob.  12.  If  the  chance  is  even  for  head  and  tail  in  tossing  a 
coin,  what  is  the  probability  that  in  a  million  throws  the  difference 
between  heads  and  tails  will  exceed  1500  ? 

ART.  7.     INVERSE  PROBABILITIES.* 

If  an  observed  event  can  be  attributed  to  any  one  of  several 
causes,  what  is  the  probability  that  any  particular  one  of  these 
causes  produced  the  event  ?  To  put  the  question  in  a  concrete 
form,  suppose  a  white  ball  has  been  drawn  from  one  of  two 
urns,  £/i  containing  3  white*  and  5  black  balls,  and  U+  contain- 
ing 2  white  and  4  black  balls  ;  and  that  the  probability  in  favor 
of  each  urn  is  required.  If  £/,  is  as  likely  to  have  been 
chosen  as  £/,,  the  probability  that  £/,  was  chosen  is  1/2.  After 
such  choice  the  probability  of  drawing  a  white  ball  from  Ul  is 
3/8.  Before  drawing,  therefore,  the  probability  of  getting  a 
white  ball  from  Ul  was  1/2  X  3/8  =  3/16,  by  Art.  5.  Similarly, 
before  drawing  the  probability  of  getting  a  white  ball  from  £/"a 
was  1/2  X  2/6  =  1/6.  These  probabilities  will  remain  un- 
changed if  the  number  of  balls  in  either  urn  be  increased  or 

*  See  Poisson,  Probabilite  des  Jugements,  pp.  81-83. 
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diminished  so  long  as  the  ratio  of  white  to  black  balls  is  kept 
constant.  Make  these  numbers  the  same  for  the  two  urns. 
Thus  let  the  first  contain  9  white  and  15  black,  and  the  second 
8  white  and  16  black ;  whence  the  above  probabilities  may  be 
written  1/2  x  9/24  and  1/2  X  8/24.  It  is  now  seen  that  there 
are  (9  -\-  8)  cases  favorable  to  the  production  of  a  white  ball, 
each  of  which  has  the  same  antecedent  probability,  namely,  1/2. 
Since  the  fact  that  a  white  ball  was  drawn  excludes  considera- 
tion of  the  black  balls,  the  probability  that  the  white  ball  came 
from  £/,  is  9/17  and  that  it  came  from  U^  is  8/17;  and  the  sum 
of  these  is  unity,  as  it  should  be. 

To  generalize  this  result,  let  there  be  m  causes,  Clf  Ct,  .  . .  Cm. 
Denote  their  direct  probabilities  by  q^  q^...qm\  their  antecedent' 
probabilities  by  rl9  ra,  . .  .  rm  ;  and  their  resultant  probabilities 
on    the   supposition   of   separate    existence    by  /,,/3,  .  .  . pm. 

That  is, 

A  =  9Si,    A  =  QSv  •  •  •  A.  =  qm*m.  (0 

Let  D  be  the  common  denominator  of  the  right-hand  mem- 
bers in  (i),  and  denote  the  corresponding  numerators  of  the 
several  fractions  by  slt  s^,  .  .  .  sm.  Then 

A  =  sJD,    A  =  sJD,  .../,  =  sJD ; 

and  it  is  seen  that  there  are  in  all  (s}  +  s*  +  •  •  •  -O  equally 
possible  cases,  and  that  of  these  ^  are  favorable  to  C19  s^'  to 
C9,  .  .  .  Hence,  if  Plf  P9,  .  .  .  Pm  denote  the  probabilities  of 
the  several  causes  on,  the  supposition  of  their  coexistence, 

Pt  =  *,/(*,  +*,+•••*,)  =  A/(A  +A  +  •  •  •  A.)- 

Thus  in  general 

P,=pj2p,     P,  =  p,/7.p,...Pm=pm/2p.  (2) 

To  illustrate  the  meaning  of  these  formulas  by  the  above 
concrete  case  of  the  urns  it  suffices  to  observe  that 
for  U19      q,  =  3/8     and     r,  =  1/2, 
for  U^      ^  =  1/3     and     r,  =  1/2  ; 
whence  /,  =  3/16,    /,  =  1/6,    /t  +  A  =  J7/48  ; 

and  Pi  =9/17,    ^  =  8/17. 

As  a  second  illustration,  suppose  it  is  known  that  a  white 
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ball  has  been  drawn  from  an  urn  which  originally  contained  nt 
balls,  some  of  them  being  black,  if  all  are  not  white.  What  is 
the  probability  that  the  urn  contained  exactly  n  white  balls? 
The  facts  are  consistent  with  m  different  and  equally  probable 
hypotheses  (or  causes),  namely,  that  there  were  I  white  and 
(;;/  —  i)  black  balls,  2  white  and  (m  —  2)  black  balls,  etc. 
Hence  in  (i),  ql  =  q^  =  .  .  .  =  I,  and 

/>,  =  i/m,    /,  =  2/m,  .  .  .  pn  =  n/m,  .  .  .  pm  -  m/m. 
Thus  2/=  (i/2)(m+  i), 

and  P= 


This  shows,  as  it  evidently  should,  that  n  =  m  is  the  most 
probable  number  of  white  balls  in  the  urn.  The  probability 
for  this  number  is  Pm  =  2/(m  -j-  i)>  which  reduces,  as  it  ought, 
to  I  for  m  —  i  . 

Formulas  (i)  and  (2)  may  also  be  applied  to  the  problem  of 
estimating  the  probability  of  the  occurrence  of  an  event  from 
the  concurrent  testimony  of  several  witnesses,  Xlf  X^  .  .  . 
Denote  the  probabilities  that  the  witnesses  tell  the  truth  by 
xl  %„  .  .  .  Then,  supposing  them  to  testify  independently,, 
the  probability  that  they  will  concur  in  the  truth  concerning 
the  event  is  xjc^  .  .  .  ;  while  the  probability  that  they  will  con- 
cur in  the  only  other  alternative,  falsehood,  is  (i  —  x^(\  —  ;ra)  .  .  . 
The  two  alternatives  are  equally  possible.  Hence  by  equations 
(i)  and  (2) 

p,  =  *,*,  .  .  .,  /„  =  (I  -  *,)(!  -*J..  ., 

•**  •  •  • 


p   = 


_ 
..  +  (I  -*,)(!  -*.)...' 

(I  -*.)(!.-*.)... 


PI  being  the  probability  for  and  P^  that  against  the  event 

To  illustrate  (3),  if  the  chances  are  3  to  i  that  X,  tells  the 
truth  and  5  to  I  that  X,  tells  the  truth,  x,  =  3/4,  x^  —  5/6,  and 
Pl=  15/16;  or,  the  chances  are  15  to  i  that  an  event  occurred 
if  they  agree  in  asserting  that  it  did.* 

*  For  some  interesting  applications  of  equations  (')  ^ee  note  E  of  Appendix 
(0  the  Ninth  Bridgewater  Treatise  by  Charles  Babbage  (London,  1838). 
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It  is  of  theoretical  interest  to  observe  that  if  Xl,  ;ra,  .  .  .  in 
(3)  are  each  greater  than  1/2,  Pl  approaches  unity  as  the 
number  of  witnesses  is  indefinitely  increased. 

Prob.  13.  The  groups  of  numbers  of  one  figure  each,  two  figures 
•each,  three  figures  each,  etc.,  which  it  is  possible  to  form  from  the 
nine  digits  i,  2,  ...  9  are  printed  on  cards  and  placed  severally  in 
nine  similar  urns.  What  is  the  probability  that  the  number  777  will 
be  drawn  in  a  single  trial  by  a  person  unaware  of  the  contents  of 
the  urns  ? 

Prob.  14.  How  many  witnesses  whose  credibilities  are  each  3/4 
are  essential  to  make  Pl  =  0.999  in  equation  (3)  ? 

ART.  8.    PROBABILITIES  OF  FUTURE  EVENTS. 

Equations  (2)  of  Art.  7  may  be  written  in  the  following 
manner: 

LJ  — _?_          LJ!  —  _i_  /T\ 

A      A          '  Pm      2p' 

If  plt  /^  ...  pm  are  found  by  observation,  Plt  Pat  .  .  .  Pm  will  ex- 
press the  probabilities  of  the  corresponding  causes  or  their 
effects.  When,  as  in  the  case  of  most  physical  facts,  the  num- 
ber of  causes  and  events  is  indefinitely  great,  the  value  of  any 
p  or  P  in  (i)  becomes  indefinitely  small,  and  the  value  of  2p 
must  be  expressed  by  means  of  a  definite  integral.  Let  x  de- 
note the  probability  of  any  particular  cause,  or  of  the  event  to 
which  it  gives  rise.  Then,  supposing  this  and  all  the  other 
causes  mutually  exclusive,  (i  —  x)  will  be  the  probability 
against  the  event.  Now  suppose  it  has  been  observed  that  in 
(m  _j_  ffy  cases  the  event  in  question  has  occurred  m  times  and 
failed  n  times.  The  probability  of  such  a  concurrence  is,  by 
Art.  5,  cxm(i  —  x)n,  where  c  is  a  constant.  Since  x  is  unknown, 
it  may  be  assumed  to  have  any  value  within  the  limits  o  and  I ; 
and  all  such  values  are  a  priori  equally  possible.  Put 

y  —  cxm(i  —  x)n- 

Then  evidently  the  probability  that  x  will  fall  within  any  as- 
signed possible  limits  a  and  b  is  expressed  by  the  fraction 


O  i 

fyte/fydz; 
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so  that  the  probability  of  any  particular  x  is  given  by 

p= 


This  may  be  regarded  as  the  antecedent  probability  of  the 
cause  or  event  in  question. 

What  then  is  the  probability  that  in  the  next  (r  -f-  s)  trials 
the  event  will  occur  r  times  and  fail  s  times,  if  no  regard  is  had 
of  the  order  of  occurrence?  If  x  were  known,  the  answer 
would  be  by  Arts.  2  and  5 


But  since  x  is  restricted  only  by  the  condition  (2),  the  required 
probability  will  be  found  by  taking  the  product  of  (2)  and  (3) 
and  integrating  throughout  the  range  of  x.  Thus,  calling  the 
required  probability  Qt 


/  xm(i  — 


The  definite  integrals  which  appear  here  are  known  as  Gamma 
functions.  They  are  discussed  in  all  of  the  higher  treatises  on 
the  Integral  Calculus.  Applying  the  rules  derived  in  such 
treatises  there  results  * 

(r  +  s)\(m  +  r)\(n  +  s)\(m  +  n+i)l 
rislml*l(m  +  n  +  r  +  s  +i)\      * 

If  regard  is  had  to  the  order  of  occurrence  of  the  event  ; 
that  is,  if  the  probability  required  is  that  of  the  event  happen- 
ing r  times  in  succession  and  then  failing  s  times  in  succession, 

*  It  is  a  remarkable  fact  that  formula  (5)  is  true  without  restriction  as  to 
values  of  m,  n,  r.  s.  The  formula  may  be  established  by  elementary  considera- 
tions, as  was  done  by  Prevost  and  Lhuilier,  1795.  See  Todhunter's  History  of 
he  Theory  of  Probability,  pp.  453~457- 
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the  factor  (r  +  s)\/rlsl  in  (3),  (4),  (5)  must  be  replaced  by 
unity. 

To  illustrate  these  formulas,  suppose  first  that  the  event 
has  happened  m  times  and  failed  no  times.  What  is  the  prob- 
ability that  it  will  occur  at  the  next  trial  ?  In  this  case  (4) 
gives 

Q  =       xm^dx         x™dx  =  (m  +  i)/(m  +  2). 


When  m  is  large  this  probability  is  nearly  unity.  Thus,  the 
sun  has  risen  without  failure  a  great  number  of  times  m  ;  the 
probability  that  it  will  rise  to-morrow  is 


which  is  practically  i. 

Secondly,  suppose  an  urn  contains  white  and  black  balls  in 
an  unknown  ratio.  If  in  ten  trials  7  white  and  3  black  balls 
are  drawn,  what  is  the  probability  that  in  the  next  five  trials 
2  white  and  3  black  balls  will  be  drawn?  The  application  of 
(5)  supposes  the  ratio  of  the  white  and  black  balls  in  the  urn 
to  remain  constant.  This  will  follow  if  the  balls  are  replaced 
after  each  drawing,  or  if  the  number  of  balls  in  the  urn  is  sup- 
posed infinite.  The  data  give 

m  =  7,  n  =  3,  r  =  2,  s  =  3, 

m  -j-  r  =  9,     n  -{-  s  =  6,     r  -f-  s  =  5,     m  -{-  n-{-  1  =  11, 
m-\-n-\-r-\-s-{-i  =  16. 

Thus  by  (5) 

5!o!6!n! 


Suppose  there  are  two  mutually  exclusive  events,  the  first 
of  which  has  happened  m  times  and  the  second  n  times  in 
m  _|_  n  trials.  What  is  the  probability  that  the  chance  of  the 
occurrence  of  the  first  exceeds  1/2  ?  The  answer  to  this  ques- 
tion is  given  directly  by  equation  (2)  by  integrating  the  nume- 
rator between  the  specified  limits  of  x.  That  is, 
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'(i  -  xfdx 


I*  xm(\  —  x)ndx 


Thus,  if  m  =  i  and  «  =  o,  P  —  3/4  ;  or  the  odds  are  three  to 
one  that  the  event  is  more  likely  to  happen  than  not.  Simi- 
larly, if  the  event  has  occurred  m  times  in  succession, 

/>=!--  (1/2)"+', 
which  approaches  unity  rapidly  with  increase  of  n. 

ART.  9.    THEORY  OF  ERRORS. 

The  theory  of  errors  may  be  defined  as  that  branch  of  math- 
ematics which  is  concerned,  first,  with  the  expression  of  the  re- 
sultant effect  of  one  or  more  sources  of  error  to  which  com- 
puted and  observed  quantities  are  subject  ;  and,  secondly,  with 
the  determination  of  the  relation  between  the  magnitude  of 
an  error  and  the  probability  of  its  occurrence.  In  the  case  of 
computed  quantities  which  depend  on  numerical  data,  such  as 
tables  of  logarithms,  trigonometric  functions,  etc.,  it  is  usually 
possible  to  ascertain  the  actual  values  of  the  resultant  errors. 
In  the  case  of  observed  quantities,  on  the  other  hand,  it  is  not 
generally  possible  to  evaluate  the  resultant  actifal  error,  since 
the  actual  errors  of  observation  are  usually  unknown.  In  either 
case,  however,  it  is  always  possible  to  write  down  a  symbolical 
expression  which  will  show  how  different  sources  of  error  enter 
and  affect  the  aggregate  error  ;  and  the  statement  of  such  an 
expression  is  of  fundamental  importance  in  the  theory  of  errors. 

To  fix  the  ideas,  suppose  a  quantity  Q  to  be  a  function  of 
several  independent  quantities  x,  y,  z  .  .  .  ;  that  is, 

Q  =/**,**...). 

and  let  it  be  required  to  determine  the  error  in  Q  due  to  errors 
in  x,  y,  z  .  .  .  Denote  such  errors  by  AQ,  Ax,  Ay,  Az  .  .  . 
Then,  supposing  the  errors  so  small  that  their  squares,  prod- 
ucts, and  higher  powers  may  be  neglected,  Taylor's  series  gives 
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&-+g<H.S*+-  (0 

This  equation  may  be  said  to  express  the  resultant  actual  error 
of  the  function  in  terms  of  the  component  actual  errors,  since 
the  actual  value  of  AQ  is  known  when  the  actual  errors  of 
x,  y,  2  .  .  .  are  known.  It  should  be  carefully  noted  that  the 
quantities  x}  y,  z  .  ..  are  supposed  subject  to  errors  which  are 
independent  of  one  another.  The  discovery  of  the  independent 
sources  of  error  is  sometimes  a  matter  of  difficulty,  and  in  general 
requires  close  attention  on  the  part  of  the  student  if  he  would 
avoid  blunders  and  misconceptions.  Every  investigator  in  work 
of  precision  should  have  a  clear  notion  of  the  error-equation  of 
the  type  (i)  appertaining  to  his  work ;  for  it  is  thus  only  that 
he  can  distinguish  between  the  important  and  unimportant 
sources  of  error. 

Prob.  15.  Write  out  the  error-equation  in  accordance  with  (i) 
for  the  function  Q  =  xyz  -f-  x*  log  (y/z). 

Prob.  1 6.  In  a  plane  triangle  a/b  —  sin  A  /sin  B.  Find  the  error 
in  a  due  to  errors  in  b,  A,  and  B. 

Prob.  17.   Suppose  in  place  of  the  data  of  problem   16  that  the 
angles  used  in  computation  are  given  by  the  following  equations  : 
^  =  ^+1(180°-^-  B,-  Q,  B  =  B,  +  t(i8o°  -  A,-B^  -  C,), 
where  Aiy  Bl ,  Cl  are  observed  values.     What  then  is  Aa* 

Prob.  1 8.  If  w  denote  the  weight  of  a  body  and  r  the  radius  of 
the  earth,  show  that  for  small  changes  in  altitude,  Aw/w=  —  2Ar/r\ 
whence,  if  a  precision  of  one  part  in  500000000  is  attainable  in  com- 
paring two  nearly  equal  masses,  the  effect  of  a  difference  in  altitude 
of  one  centimeter  in  the  scale-pans  of  a  balance  will  be  noticeable.* 

ART.  10.     LAWS  OF  ERROR. 

A  law  of  error  is  a  function  which  expresses  the  relative 
frequency  of  occurrence  of  errors  in  terms  of  their  magnitudes. 
Thus,  using  the  customary  notation,  let  e  denote  the  magni- 

*  This  problem  arose  with  the  International  Bureau  of  Weights  and  Measures, 
whose  work  of  intercomparison  of  the  Prototype  Kilogrammes  attained  a  pre- 
cision indicated  by  a  probable  error  of  1/500  ooo  oooth  part  of  a  kilogramme. 
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tude  of  any  error  in  a  system  of  possible  errors.  Then  the  law 
of  such  system  may  be  expressed  by  an  equation  of  the  form 

y  =  0(<0-  (0 

Representing  e  as  abscissa  and  y  as  ordinate,  this  equation 
gives  a  curve  called  the  curve  of  frequency,  the  nature  of  which, 
as  is  evident,  depends  on  the  form  of  the  function  0.  This 
equation  gives  the  relative  frequency  of  occurrence  of  errors  in 
the  system  ;  so  that  if  e  is  continuous  the  probability  of  the 
occurrence  of  any  particular  error  is  expressed  byyde  =  <f>(e)de\ 
which  is  infinitesimal,  as  it  plainly  should  be,  since  in  any  con- 
tinuous system  the  number  of  different  values  of  e  is  infinite. 

Consider  the  simplest  form  of  0(e),  namely,  that  in  which 
0(e)  =  c,  a  constant.  This  form  of  0(e)  obtains  in  the  case  of 
the  errors  of  tabular  logarithms,  natural  trigonometric  func- 
tions, etc.  In  this  case  all  errors  between  minus  a  half-unit 
and  plus  a  half-unit  of  the  last  tabular  place  are  equally  likely 
to  occur.  Suppose,  to  cover  the  class  of  cases  to  which  that 
just  cited  belongs,  all  errors  between  the  limits  —  a  and  -\-a 
are  equally  likely  to  occur.  The  probability  of  any  individual 
error  will  then  be  <f)(e)de  —  cde,  and  the  sum  of  all  such  prob- 
abilities, by  equation  (2),  Art.  4,  must  be  unity.  That  is, 


«  « 

J*(f>(e)de  =  cCde  =  I. 


(2) 


This  gives  c  —  1/20,  or  by  (i)  y  =  i/2a.     The  curve  of  fre- 
quency in  this  case  is  shown  in  the  figure, 
AB  being  the  axis  of  e  and  OQ  that  of/. 
It  is  evident  from  this  diagram  that  if  the 
errors   of  the  system  be  considered  with 
respect  to  magnitude   only,  half  of  them 
should  be  greater  and  half  less  than  a/2. 
'  This  is  easily  found  to  be  so  in  the  case  of 
tabular  logarithms,  etc. 

As  a  second  illustration  of  (i),  suppose  y  and  e  connected 
by  the  relation^  =  c  Va*  —  e*,  where  a  is  the  radius  of  a  circle, 
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c  a  constant,  and  e  may  have  any  value  between  —  a  and  -f  «. 
Then  the  condition 

+  a 

cj'de  Va^7'  =  i 

—  a 

gives  c  =  2/(a*n).  In  this,  as  in  the  preceding  case,  0(+  e)  = 
0(—  e),  the  meaning  of  which  is  that  positive  and  negative 
errors  of  the  same  magnitude  are  equally  likely  to  occur.  It 
will  be  noticed,  however,  that  in  the  latter  case  small  errors 
have  a  much  higher  probability  than  those  near  the  limit  a+ 
while  in  the  former  case  all  errors  have  the  same  probability. 

In  general,  when  e  is  continuous  0(e)  must  satisfy  the  condi- 
tion J  (p(e)de  =  I,  the  limits  being  such  as  to  cover  the  entire 

range  of  values  of  e.  The  cases  most  commonly  met  with  are 
those  in  which  0(e)  is  an  even  function,  or  those  in  which 
0(+  e)  =  0(—  e).  In  such  cases,  if  ±  a  denote  the  limiting 
value  of  e, 

-fa  a 

eXe  =  i.  (3) 


ART.  11.    TYPICAL  ERRORS  OF  A  SYSTEM. 

Certain  typical  errors  of  a  system  have  received  special 
designations  and  are  of  constant  use  in  the  literature  of  the 
theory  of  errors.  These  special  errors  are  the  probable  error, 
the  mean  error,  and  the  average  error.  The  first  is  that  error 
of  the  system  of  errors  which  is  as  likely'to  be  exceeded  as 
not  ;  the  second  is  the  square  root  of  the  mean  of  the  squares 
of  all  the  errors  ;  and  the  third  is  the  mean  of  all  the  errors 
regardless  of  their  signs.  Confining  attention  to  systems  in 
which  positive  and  negative  errors  of  the  same  magnitude 
are  equally  probable,  these  typical  errors  are  defined  mathe- 
matically as  follows.  Let 

ep  =  the  probable  error, 

em=  the  mean  error, 

ea  =  the  average  error. 
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Then,  observing  (2),  of  Art.  10, 
f<t>(e)de  =f(t>(e)de  = 


e   = 


(0 


The  student  should  seek  to  avoid  the  very  common  misap- 
prehension of  the  meaning  of  the  probable  error.  It  is  not 
*'  the  most  probable  error,"  nor  "  the  most  probable  value  of 
the  actual  error";  but  it  is  that  error  which,  disregarding  signs, 
•would  occupy  the  middle  place  if  all  the  errors  of  the  system 
-were  arranged  in  order  of  magnitude.  A  few  illustrations  will 
suffice  to  fix  the  ideas  as  to  the  typical  errors.  Thus,  take  the 
simple  case  wherein  <j)(e)  =  c=  1/2^,  which  applies  to  tabular 
logarithms,  etc.  Equations  (i)  give  at  once 

€f=  ±  -a,     em  =  ±  -  ^3~,     ea  =  ±  -a. 

For  the  case  of  tabular  values,  a  —  0.5  in  units  of  the  last 
"tabular  place.  Hence  for  such  values 

e,  =  ±0.25,     ew=±o.29,     €,,=  ±0.25. 

Prob.  19.  Find  the  typical  errors  for  the  cases  in  which  the  law 
of  error  is  0(e)  =  cVa*  —  e\  0(e)  =  d±a^e),  (f)(e)=c  cosa(7re  /2<z); 
£  being  a  constant  to  be  determined  in  each  case  and  e  having  any 
value  between  —  a  and  +  &• 

ART.  12.    LAWS  OF  RESULTANT  ERROR. 

When  several  independent  sources  of  error  conspire  to  pro- 
duce a  resultant  error,  as  specified  by  equation  (i)  of  Art.  9, 
there  is  presented  the  problem  of  determining  the  law  of  the 
resultant  error  by  means  of  the  laws  of  the  component  errors. 
The  algebraic  statement  of  this  problem  is  obtained  as  follows 
for  the  case  of  continuous  errors  : 

In  the  equation  (i),  Art.  9,  write  for  brevity 
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and  let  the  laws  of  error  of  e,  elt  ea,  .  .  .  be  denoted  by  0(e)r 
Pi(ei)>  0a(ea)  •  •  •     Then  the  value  of  e  is  given  by 

*=*,  +  *,.+•••  (i> 

The  probabilities  of  the  occurrence  of  any  particular  values- 
of  el9  ea,  .  .  .  are  given  by  0,(eiXei»  '  0a(eaKea»  •  •  •  ;  and  the 
probability  of  their  concurrence  is  the  probability  of  the  cor- 
responding value  of  e.  But  since  this  value  may  arise  in  an 
infinite  number  of  ways  through  the  variations  of  e,,  ea,  .  .  » 
over  their  ranges,  the  probability  of  e-,  or  </>(e)de,  will  be 
expressed  by  the  integral  of  01(e1)^e102(eaye2  .  .  .  subject  to- 
the  restriction  (i).  This  latter  gives  e,  =  e—  ea  —  e,  .  .  .,  and 
de^  =  de  for  the  multiple  integration  with  respect  to  ea,  e,,  .  .  „ 
Hence  there  results 


or 


It  is  readily  seen  that  this  formula  will  increase  rapidly  in 
complexity  with  the  number  of  independent  sources  of  error.*" 
For  some  of  the  most  important  practical  applications,  how- 
ever, it  suffices  to  limit  equation  (2)  to  the  case  of  two  inde- 
pendent sources  of  error,  each  of  constant  probability  within 
assigned  limits.  Thus,  to  consider  this  case,  let  e1  vary  over 
the  range  —  a  to  -f-  a,  and  e2  vary  over  the  range  —  b  to  +  b* 
Then  by  equation  (2),  Art.  10, 

0,^)  -  i/(2*)>     0,(O  = 
Hence  equation  (2)  becomes 


In    evaluating    this   integral    ea   must    not    surpass    ±  b   and' 

6j  =  e  —  e2  must  not  surpass  ±  a.     Assuming  a  >  b,  the  limits 

of   the  integral  for  any  value  of   €=€,  +  €,    lying   between 

-  (a  +  b)  and  —  (a  —  b)  are  —  b  and  +  (e  +  a).     This  fact  is- 

*  The  reader  desirous  of  pursuing  this  phase  of  the  subject  should  consult 
Bessel's  Untersuchungen  ueber  die  Wahrscheinlichkeit  der  Beobachtungsfehler; 
Abhandlungen  von  Bessel  (Leipzig,  1876),  Vol.  II. 
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made  plain  by  a  numerical  example.  For  instance,  suppose 
a  =  5  and  b  =  3.  Then  —  (a  +  b)  =  —  8  and  —  (a  —  b)  —  —  3. 
Take  e  =  —  6,  a  number  intermediate  to  —  8  and  —  3.  Then 
the  following  are  the  possible  integer  values  of  e,  and  ea  which 
\vill  produce  e  =  —  6 : 

e  e,      e,  limits  of  e, 

-6=  -5-  i,       -  i  =  +  (e  +  a)t 
=  -4-2, 

=  -  3  -  3,     -  3  =  -  *• 

Similarly,  the  limits  of  ea  for  values  of  e  lying  between 
—  (a  —  b)  and  -\-  (a  —  b)  are  —  b  and  -\-  b ;  and  the  limits  of 
€9  for  values  of  e  between  -\-  (a  —  b)  and  -|-  (a  -f-  b)  are  -}-  (e  —  a) 
and  -f-  b.  Hence 


-  (3) 


Thus  it  appears  that  in  this  case  the  graph  of  the  resultant 

law  of  error  is  represented  by  the  upper  base  and  the  two  sides 

of   a  trapezoid,  the  lower  base  being  the 

axis  of  e  and  the  line  joining  the  middle 

points  of  the  bases  being  the  axis  of  0(e). 

HIIOIIII  .,          .  .          *  .        _, 

(See  the  first  figure  in  Art.  13.)     The  prop- 
in  1 101 1  ill  . 

erties  of  (3),  including  the  determination 

of  the  limits,  are   also  illustrated   by  the 
IIIIIIIOIIIIIII  .,       . 

adjacent  trapezoid  of    numerals   arranged 

IIIIIIIIOIIIIIIII         J 

to  represent  the  case  wherein  a  =  0.5  and 

£  =  0.3.     The  vertical  scale,  or  that  for  0(e),  does  not,  how- 
ever, conform  exactly  to  that  for  e. 


ERRORS  OF  INTERPOLATED  VALUES.  37 

Prob.  20.  Prove  that  the  values  of  0(e)  as  given  by  equation  (3) 
satisfy  the  condition  specified  in  equation  (3),  Art.  10. 

Prob.  21.  Examine  equations  (3)  for  the  cases  wherein  a  =  £and 
b  =  o;  and  interpret  for  the  latter  case  the  first  and  last  of  (3). 

Prob.  22.  Find  from  (3),  and  (i)  of  Art.  n,  the  probable  error  of 
the  sum  of  two  tabular  logarithms. 

ART.  13.  ERRORS  OF  INTERPOLATED  VALUES. 
Case  I. — One  of  the  most  instructive^cases  to  which  formulas 
(3)  of  Art.  12  are  applicable  is  that  of  interpolated  logarithms, 
trigonometric  functions,  etc.,  dependent  on  first  differences. 
Thus,  suppose  that  vt  and  v^  are  two  tabular  logarithms,  and 
that  it  is  required  to  get  a  value  v  lying  /  tenths  of  the  interval 
from  vl  towards  vt.  Evidently 

v  =  v*  +  fa  -  «0  t=(i  —  t}v,  +  tv%  ; 

and  hence  if  e,  elt  e^  denote  the  actual  errors  of  v,  vltvt  ,  re- 
spectively, 

*  =  0  -')'.  +  ",-  (I) 

It  is  to  be  carefully  noted  here  that  e  as  given  by  (i)  re- 
quires the  retention  in  v  of  at  least  one  decimal  place  be- 
yond the  last  tabular  place.  For  example,  let  v  =  log  (24373) 
from  a  5-place  table.  Then  v^  =  4.38686,  vt  =  4.38703, 
^  —  Vt=  +0.00017,  /  =  0.3,  and  v  =4.38691.1.  Likewise,  as 
found  from  a  7-place  table,  el  =  —  0.45,  et  =  -f-  0.37  in  units  of 
the  fifth  place;  and  hence  by  (i)  e  =  —0.20.  That  is,  the 
actual  error  of  v  =  4.38691.1  is  =  0.20,  and  this  is  verified  by 
reference  to  a  7-place  table. 

The  reader  is  also  cautioned  against  mistaking  the  species 
of  interpolated  values  here  considered  for  the  species  common- 
ly used  by  computers,  namely,  that  in  which  the  interpolated 
value  is  rounded  to  the  nearest  unit  of  the  last  tabular  place. 
The  latter  species  is  discussed  under  Case  II  below. 

Confining  attention  now  to  the  class  of  errors  specified  by 
equation  (i),  there  result  in  the  notation  of  the  preceding 

article 

eI  =  (i—  /X,     ea  =  tea,     and     e  =  e  =  e,  +  e2 ; 

and  since  el  and  e^  each  vary  continuously  between  the  limits 
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it  0.5  of  a  unit  of  the  last  tabular  place,  a  and  b  in  equations 
(3)  of  that  article  have  the  values 

a  —  0.5(1  -  /),     b=o.$t. 

Hence  the  law  of  error  of  the  interpolated  values  is  ex- 
pressed as  follows  : 


-     ^L_  for  values  of  e  betw.  —0.5  and  —  (0.5—  /),  1 


=  -    — for  values  of  e  betw.  —  (0.5—  t)  and  +(0.5— -/), 
i  —  * 

=  -r-        —  for  values  of  e  betw.  +(0.5—  t)  'and  +0.5. 


(2) 


The  graph  of  0(e)  for  /  =  1/3  is  shown  by  the  trapezoid 
AB,  BC,  CD  in  the  figure  on  page  40.  Evidently  the  equa- 
tions (2)  are  in  general  represented  by  a  trapezoid,  which  degen- 
erates to  an  isosceles  triangle  when  /  =  1/2. 

The  probable,  mean,  and  average  errors  of  an  interpolated 
value  of  the  kind  in  question  are  readily  found  from  (2),  and 
from  equations  (i)  of  Art.  II,  to  be 

-/  for      o</<  i 


=  1/2  -  (i/2)V2t(i  -  /)      for  1/3  <  /  <  2/3, 
=  i/4/  for  2/3  </<  i. 


e«  = 


96(1  _/)/ 


(3) 


- 

gr          ^  1/2  </<i.. 

It  is  thus  seen  that  the  probable  error  of  the  interpolated 
value  here  considered  decreases  from  0.25  to  0.15  of  a  unit  of 
the  last  tabular  place  as  t  increases  from  o  to  0.5.  Hence  such 
values  are  more  precise  than  tabular  values  ;  and  the  computer 
who  desires  to  secure  the  highest  attainable  precision  with  a 
given  table  of  logarithms  should  retain  one  additional  figure 
beyond  the  last  tabular  place  in  interpolated  values. 
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Case  II.  —  Recurring  to  the  equation  v  =  vv  +  t(v^  —  v>)  for  an 
interpolated  value  v  in  terms  of  two  consecutive  tabular  values 
vl  and  ?>a,  it  will  be  observed  that  if  the  quantity  t(vt  —  z/,)  is 
rounded  to  the  nearest  unit  of  the  last  tabular  place,  a  new  error 
is  introduced.  For  example,  if  ^'=log  1633  =  3.21299,  and 
v,  =  log  1634  =  3.21325  from  a  5-place  table,  v,  -  vi  =  +  26 
units  of  the  last  tabular  place  ;  and  if  t  =  1/3,  t(v9  —  v^  =  8f  ; 
so  that  by  the  method  of  interpolation  in  question  there  results 
v  =  3.21299  +  9  =  3.21308.  Now  the  actual  errors  of  vl  and 
£>8  are,  as  found  from  a  /-place  table,  —  0.38  and  +0.21  in  units 
of  the  fifth  place.  Hence  the  actual  error  of  v  is  by  equation 
0)»  f  X  -  0.38  +  £  X  +  0.21  —  J  =  -  0.52,  as  is  shown  di- 
rectly by  a  7-place  table. 

It  appears,  then,  that  in  this  case  the  error-equation  cor- 
responding to  (i)  is 

*=  (I  -*X  +  ",+'.,  (4) 

wherein  el  and  e^  are  the  same  as  in  (i)  and  ^3is  the  actual  error 
that  comes  from  rounding  t(vt  —  #,)  to  the  nearest  unit  of  the 
last  tabular  place. 

The  error  ez,  however,  differs  radically  in  kind  from  ^,  and 
<?3.  The  two  latter  are  continuous,  that  is,  they  may  each  have 
any  value,  between  the  limits  —  0.5  and  +0.5  ;  while  ea  is  dis- 
continuous, being  limited  to  a  finite  number  of  values  depend- 
ent on  the  interpolating  factor  /.  Thus,  for  /  =  1/2  the  only 
possible  values  of  ea  are  o  -f-  1/2,  and  —  1/2  ;  likewise  for  /  = 
1/3,  the  only  possible  values  of  es  are  O,  -j-  1/3,  and  —  1/3.  It 
is  also  clear  that  the  maximum  value  of  ^,  which  is  constant  and 
equal  to  1/2  for  (i),  is  variable  for  (4)  in  a  manner  dependent 
on  /.  For  example,  in  (4), 

The  maximum  of  e  =  1/2  +  1/2  =  i,  for  t  =  1/2, 
«  ,  =  1/2  +  1/3  =  5/6,  "  /  =  1/3, 
"  f  =  1/2  +  1/2  =  i,  "  /  =  1/4, 

"  "  e  =   1/2  +  2/5  =  9/10  "      /  =:   1/5. 

The  determination  of  the  law  of  error  for  this  case  presents 
some  novelty,  since  it  is  essential  to  combine  the  continuous 
errors  (i  —  t)e^  and  te^  with  the  discontinuous  error  ey  The 


6      B 
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simplest  mode  of  attacking  the  problem  seems  to  be  the  fol- 
lowing quasi-geometrical  one.     In  the  notation  of  Arts.  12  and 
13,  put  in  (4)  e  —  e,  (i  —  t)el  =  e,,  /^a  =  ea,  and  e,  =  e3.     Then 
e  =  (e,  -}-  e,)  +  e,.  (5) 

The  law  of  error  for  (el  -f-  ea)  is  given  by  equation  (2)  for  any 
value  of  /.  Hence  for  a  given  value  of  /  there  will  be  as  many 
expressions  of  0(e)  as  there  are  different  values  of  es.  The 
graphs  of  0(e)  will  all  be  of  the  same  form  but  will  be  differently 
placed  with  reference  to  the  axis  of  0(e).  Thus,  if  /  =  1/3  the 

values  of  e8  are  —  1/3,  o,  and 
-f-  i/3,  and  these'  are  equally 
likely  to  occur.  For  e8  =  o  the 
graph  is  given  directly  by  (2), 
and  is  the  trapezoid  ABCD 
symmetrical  with  respect  to  OQ. 
For  e3  =  —  1/3  the  graph  is 
abQd,  of  the  same  form  as 
ABCD  but  shifted  to  the  left 
by  the  amount  of  es  =  —  1/3. 
Similarly,  the  graph  for  the  case 
of  e,  =  +  1/3  is  a'Qb'd',  and  is  produced  by  shifting  ABCD  to 
the  right  by  an  amount  equal  to  +  1/3. 

Now,  since  the  three  systems  of  errors  for  this  case  are 
equally  likely  to  occur,  they  may  be  combined  into  one  system 
by  simple  addition  of  the  corresponding  element  areas  of  the 
several  graphs.  Inspection  of  the  diagram  shows*  that  the 
resultant  law  of  error  is  expressed  by 

0(e)  =  (i/4)(5  +  6e)     for  -  5/6  <  e  <  -  1/6,  <| 

=  i  for  -  1/6  <  e  <  +  1/6,  I        (6) 

=  (i/4)(5  -  6e)      for  +  1/6  <  e  <  +  5/6.  J 

This  is  represented  by  a  trapezoid  whose  lower  base  is  10/6, 

upper  base  2/6,  and  altitude  i. 


*  Sum  the  three  areas  and    divide   by  3  to  make   resultant   area  =  i,    as 
required  by  equation  (3),  Art.  10. 
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As  a  second  illustration,  consider  equation  (5)  for  the  case 
t  =  1/2.     In  this  case  e,  must  be  either  o  or  1/2,  the  sign  of 
which  latter  is  arbitrary.     For  e3  =  o,  equations  (2)  give 
0(e)  =  2  +  4e     for     _  1/2  <  6  <  o,  i 

=  2—  4€     for     Q-.<e<+l/2.  )  (7) 

This   function   is  represented  by  the  isosceles  triangle  A  QE 

whose  altitude  OQ  is  twice  the  base  AE. 

Similarly  0(e)  for  e,  =  -f-  1/2  would 

be  represented  by  the  triangle  AQE  dis- 
placed to  the  right  a  distance  1/2  ;  and 
if  the  two  systems  for  ea  =  o  and  ea  = 
-f-  1/2  be  combined  into  one  system, 
their  resultant  law  of  error  is  evidently 
<t>(e)  =  i+2e  for  — 1/2 «?<  o,  \ 

=  I  for  o<  e<  +  i/2,  (8) 

=  2  —  2e  for  +  1/2  <  e<  i  ;  ) 


\ 


\ 


the  graph  of  which  is  ABCD.    On  the 

other  hand,  if  the  errors  in  this  combined  system  be  considered 

with  respect  to  magnitude  only,  the  law  of  error  is 

0(e)  =  2(1  -  e)     for     o  <  e<  I,  (9) 

the  graph  of  which  is  OQD. 

The  student  should  observe  that  (6),  (7),  (8),  and  (9)  satisfy 
the  condition  J  cf>(e)de  =  i  if  the  integration  embraces  the 
whole  range  of  e. 

The  determination  of  the  general  form  of  0(e)  in  terms  of 
the  interpolating  factor  /  for  the  present  case  presents  some 
difficulties,  and  there  does  not  appear  to  be  any  published  solu- 
tion of  this  problem.*  The  results  arising  from  one  phase  of 
the  problem  have  been  given,  however,  by  the  author  in  the 
Annals  of  Mathematics,f  and  may  be  here  stated  without 
proof.  The  phase  in  question  is  that  wherein  /'is  of  the  form 
i/n,  n  being  any  positive  integer  less  than  twice  the  greatest 

*  The  author  explained  a  general  method  of  solution  in  a  paper  read  at  the 
summer  meeting  of  the  American  Mathematical  Society,  August,  1895. 
f  Vol.  II,  pp.  54-59. 
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tabular  difference  of  the  table  to  which  the  formulas  are  ap- 
plied. For  this  restricted  form  of  t  the  possible  maximum 
value  of  e  as  given  by  equation  (5)  is,  in  units  of  the  last 
tabular  place,  (2n  —  i)/n  for  n  odd  and  I  for  n  even. 

The  possible  values  of  e3  of  equation  (5)  are 

I  2  n  —  i 

O,     ±  -,     ±  -,  .  .  .  ±  —^-  for  n  odd, 

I  2  n  —  2  i 

°'    ±  »•    *?•"•*  1JT-    ±2    for  "  even- 

An  important  fact  with  regard  to  the  error  1/2  for  n  even 
is  that  its  sign  is  arbitrary,  or  is  not  fixed  by  the  computation 
as  is  the  case  with  all  the  other  errors.  However,  the  com- 
puter's rule,  which  makes  the  rounded  last  figure  of  an  inter- 
polated value  even  when  half  a  unit  is  to  be  disposed  of,  will, 
in  the  long-run,  make  this  error  as  often  plus  as  minus. 

The  laws  of  error  which  result  are  then  as  follows  : 

For  n  odd. 
0(e)  —  i  for  e  between  —  \/2n  and  ~\-  i/2n, 

0(e)  =  -  -  ±  ej  for  ebetw.  =F  i/2n  and  ^f(2n—  i)/2». 

For  n  even. 

0(e)  =  —  -  A—        -  ±  e  )  for  e  between  o  and  =F  i/«, 


e\  for  e  betw.  =f  \/n  and  =F  (»  — 


—  r(i  ±  e)          for  e  between  qp  («—  i/«)  and  q=  i. 

By  means  of  these  formulas  and  (i)  of  Art.  u  the  probable, 
mean,  and  average  errors  for  any  value  of  n  can  be  readily 
found.  The  following  table  contains  the  results  of  such  a  com- 
putation for  values  of  n  ranging  from  i  to  10.  The  maximum 
actual  error  for  each  value  of  n  is  also  added.  The  verifica- 
tion of  the  tabular  quantities  will  afford  a  useful  exercise  to  the 
student. 
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TYPICAL  ERRORS  OF  INTERPOLATED  LOGARITHMS,  ETC. 


Interpolating 
Factor. 
t=.  i/n 

Probable  Error. 
€P 

Mean  Error. 
em 

Average  Error. 
ea 

Maximum 
Actual  Error. 

l 

0.250 

0.289 

0.250 

1/2 

1/2 

.292 

.408 

•333 

I 

1/3 

.256 

•  347 

.287 

5/6 

i/4 

.276 

.382 

oi3 

i 

1/5 

.268 

•  370 

•303 

9/10 

1/6 

=  277 

.385 

.315 

r 

1/7 

.274 

.380 

•311 

13/14 

1/8 

.279 

.389 

.318 

i 

1/9 

.278 

.386 

.316 

17/18 

1/10 

.281 

.392 

.320 

i 

When  the  interpolating  factor  /  has  the  more  general  form 
m/n,  wherein  m  and  n  are  integers  with  no  common  factor,  the 
possible  values  of  es  are  the  same  as  for  /  =  i/n.  But  equa- 
tions (3)  of  Art.  12  are  not  the  same  for  t  —  m/n  as  for/  =  i/n, 
and  hence  for  the  more  general  form  of  /,  0(e)  assumes  a  new 
type  which  is  somewhat  more  complex  than  that  discussed 
above.  The  limits  of  this  work  render  it  impossible  to  extend 
the  investigation  to  these  more  complex  forms  of  0(e).  It  may 
suffice,  therefore,  to  give  a  single  instance  of  such  a  function, 
namely,  that  for  which  t  =  2/5.  For  this  case 

0(e)  =i  for  e  between  o  and  =f  i/io, 

—  (5/6X13/10  ±  e)  for  e  between  =F  i/io  and 
=  (5/3X4/5  ±  e)  for  e  between  qp  3/10  and 
=  (5/6)(9/10  ±  e)  f°r  e  between 

The  graph  of  the  right-hand  half  of  A 
this  function  is  shown  in  the  accompany- 
ing  diagram,    the   whole   graph    being 
symmetrical  with  respect  to  OA,  or  the 
axis  of  0(e). 

Attention  may  be  called  to  the  strik- 
ing resemblance  of  this  graph  to  that  of 
the  law  of  error  of  least  squares. 


3/10, 
7/io, 
9/10. 


Prob.  23.  Show  from  equations  (3)  that  em  varies  from 
=  0.29  — ,  for  /  =  o,  to  1/^24  =  0.20  -f,  for  /  =  0.5  ;  and  that  ea 
varies  from  0.25  to  i/6  for  the  same  limits. 
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Prob.  24.  Show  that  the  probable,  mean,  and  average  errors 
for  the  case  of  /=  2/5  cited  above  (p.  43)  are  ±  0.261,  ±  0.251, 
and  ±  0.290,  respectively. 

ART.  14.  STATISTICAL  TEST  OF  THEORY. 
A  statistical  test  of  the  theory  developed  in  Art.  13  may 
be  readily  drawn  from  any  considerable  number  of  actual  er- 
rors of  interpolated  values  dependent  on  the  same  interpolating 
factor.  The  application  of  such  a  test,  if  carried  out  fully  by 
the  student,  will  go  far  also  towards  fixing  clear  notions  as  to 
the  meaning  of  the  critical  errors. 

Consider  first  the  case  in  which  an  interpolated  value  falls 
midway  between  two  consecutive  values,  and  suppose  this 
interpolated  value  retains  two  additional  figures  beyond  the 
last  tabular  place.  Then  by  equations  (2),  Art.  13,  the  law  of 
error  of  this  interpolated  value  is 

0(e) '=  2  -f-  4e  for  e  between  —  0.5  and  o 
=  2  —  46  for  e  between  o  and  +  0.5. 

Hence  by  equation  (i)  of  Art.  1 1,  or  equation  (3)  of  Art.  12,  the 
probable  error  in  this  system  of  errors  is  ^  —  (J)  4/2  =  0.15. 
It  follows,  therefore,  that  in  any  large  number  of  actual  errors 
of  this  system,  half  should  be  less  and  half  greater  than  0.15. 
Similarly,  of  the  whole  number  of  such  errors  the  percentage 
falling  between  the  values  o.o  and  0.2  should  be 

+  0-2  +0.8 

J  <j)(e)de  =  2J  (2  —  4e)de  =  0.64 ; 

-0.2  0 

that  is,  sixty-four  per  cent  of  the  errors  in  question  should  be 
less  numerically  than  0.2. 

To  afford  a  more  detailed  comparison  in  this  case,  the  act- 
ual errors  of  five  hundred  interpolated  values  from  a  5-place 
table  have  been  computed  by  means  of  a  /-place  table.  The 
arguments  used  were  the /following  numbers:  20005,  2OC35> 
20065,  20105,  20135,  etc.,  in  the  same  order  to  36635.  The 
actual  and  theoretical  percentages  of  the  whole  number  of 
errors  falling  between  the  limits  o.o  and  o.i,  o.i  and  0.2,  etc., 
are  shown  in  the  tabular  form  following: 
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Limits  of  Errors.  D  Actual  Theoretical 

Percentage.  Percentage. 

o.oando.1 33.2  36 

o.  i  and  0.2 30.2  28 

0.2  and  0.3 19.0  20 

0.3  and  0.4 ,....  13.2  12 

0.4  and  0.5 4.4  4 

o.oando.15 51.4  50 

The  agreement  shown  here  between  the  actual  and  theoretical 
percentages  is  quite  close,  the  maximum  discrepancy  being  2.8 
and  the  average  1.5  per  cent. 

Secondly,  consider  the  case  of  interpolated  mid-values  of  the 
species  treated  under  Case  II  of  Art.  13.  The  law  of  error  for. 
this  case  is  given  by  the  single  equation  (9)  of  Art.  13,  namely, 
0(e)  =  2(1  —  e),  no  regard  being  paid  to  the  signs  of  the  errors. 
The  probable  error  is  then  found  from 


whence ~ep  =  i  —  \  V2  =  0.29.  Similarly,  the  percentage  of 
the  whole  number  of  errors  which  may  be  expected  to  lie,  for 
example,  between  o.o  and  0.2  in  this  system  is 


V.I 

>  /  (i  —  e)de  =  0.36. 


Using  the  same  five  hundred  interpolated  values  cited 
above,  but  rounding  them  to  the  nearest  unit  of  the  last  tabu- 
lar place  and  computing  their  actual  errors  by  means  of  a  /-place 
table,  the  following  comparison  is  afforded  : 

Actual  Theoretical 

Limits  of  Errors.  Percentage.         Percentage. 

o.o  and  0.2 35-8  36 

0.2  and  0.4 27.8  28 

0.4  and  0.6 18.6  20 

0.6ando.8 ,..  12.2  12 

0.8  and  i.o 5-6  4 

o.o  and  0.29 49-8  5° 
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The  agreement  shown  here  between  the  actual  and  theoretical 
percentages  is  somewhat  closer  than  in  the  preceding  case,  the 
maximum  discrepancy  being  only  1.6  and  the  average  only  0.6 
per  cent. 

Finally,  the  following  data  derived  from  one  thousand  act- 
ual errors  may  be  cited.  The  errors  of  one  hundred  inter- 
polated values  rounded  to  the  nearest  unit  of  the  last  tabular 
place  were  computed  *  for  each  of  the  interpolating  factors 
o.i,  0.2,  .  .  .  0.9.  The  averages  of  these  several  groups  of  act- 
ual errors  are  given  along  with  the  corresponding  theoretical 
errors  in  the  parallel  columns  below: 

Interpolating  Actual  Theoretical 

Factor.  Average  Error.      Average  Error. 

o.i 0.338  0.320 

0.2 0.288  0.303 

0.3 0.321  0.304 

0.4 O.268  0.290 

0-5 0-324  0.333 

0.6 ...    0.276  0.290 

0.7 0.321  0.304 

0.8 0.289  0.303 

0.9 0.347  0.320 

The  average  discrepancy  between  the  actual  and  theoret- 
ical values  shown  here  is  0.017.  It  is,  perhaps,  somewhat 
smaller  than  should  be  expected,  since  the  computation  of  the 
actual  errors  to  three  places  of  decimals  is  hardly  warranted 
by  the  assumption  of  dependence  on  first  differences  only. 

The  average  of  the  whole  number  of  actual  errors  in  this 
case  is  0.308,  which  agrees  to  the  same  number  of  decimals 
with  the  average  of  the  theoretical  errors,  f 

*  By  Prof.  H.  A.  Howe.  See  Annals  of  Mathematics,  Vol.  Ill,  p.  74. 
The  theoretical  averages  were  furnished  to  Prof.  Howe  by  the  author. 

f  The  reader  who  is  acquainted  with  the  elements  of  the  method  of  least 
squares  will  find  it  instructive  to  apply  that  method  to  equation  (i),  Art.  13, 
and  derive  the  probable  error  of  e.  This  is  frequently  done  without  reserve  by 
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Prob.  25.  Apply  formulas  (3)  of  Art.  12  to  the  case  of  the  sum 
or  difference  of  two  tabular  logarithms  and  derive  the  correspond- 
ing values  of  the  probable,  mean,  and  average  errors.  The  graph 
of  0(e)  is  in  this  case  an  isosceles  triangle  whose  base,  or  axis  of  e, 
is  2,  and  whose  altitude,  or  axis  of  0(e),  is  i. 

those  familiar  with  least  squares.     Thus,  the  probable  error  of  e\  or  f9  being 
0.25,  the  probable  error  of  e  is  found  to  be 


0.25     i  — 

This  varies  between  0.25  for  t  =  o  and  o.  18  for  /  =  £  ;  while  the  true  value  of 
the  probable  error,  as  shown  by  equations  (3),  Art.  13,  varies  from  0.25  to  0.15 
for  the  same  values  of  /.  It  is,  indeed,  remarkable  that  the  method  of  least 
squares,  which  admits  infinite  values  for  the  actual  errors  ei  and  e^,  should  give 
so  close  an  approximate  formula  as  the  above  for  the  probable  error  of  e. 

Similarly,  one  accustomed  to  the  method  of  least  squares  would  be  inclined 
to  apply  it  to  equation  (4),  Art.  13,  to  determine  the  probable  error  of  e.  The 
natural  blunder  in  this  case  is  to  consider  e\t  ei,  and  <?3  independent,  and  ea  like 
e\  and  e^  continuous  betweer  the  limits  o.o  and  0.5  ;  and  to  assign  a  probable 
error  of  0.25  to  each.  In  t'.is  manner  the  value 


0.25  ^2(1  -  t+t*) 

is  derived.  But  this  is  absurd,  since  it  gives  0.25  ^2  instead  of  0.25  for  /  =  o. 
The  formula  fails  then  to  give  even  approximate  results  except  for  values  of  / 
near  0.5. 
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